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Abstract 

We discuss a simple derivation of the real-time AdS / CFT prescription as an analytic continuation 
of the corresponding problem in Euclidean signature. We then extend the formalism to spinor 
operators and apply it to the examples of real-time fermionic correlators in CFTs dual to pure AdS 
and the BTZ black hole. 
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I. INTRODUCTION 



The AdS / CFT correspondence has provided important tools for studying strongly 

coupled systems by relating them to classical gravity (or string) systems. For understand- 
ing properties of such a system at finite temperature or finite density, real-time two-point 
correlation functions of composite operators are of great importance, as they encode infor- 
mation about collective behavior of the system such as transport, (the possible existence of) 
quasi-particles, etc. 

A simple prescription for calculating retarded two-point functions in AdS/CFT was first 
proposed by Son and Starinets some years ago in JJ . This prescription, which was justified 
in different ways in the literature (see e.g. 5, 6, 7, ^, {J), has been instrumental in extracting 
many important insights into strongly interacting many-body systems from AdS/CFT. 



In (lOj we observed that the prescription of [4| could be reformulated in terms of bound- 
ary values of the canonical momenta of bulk fields by treating the AdS radial direction as 
"time" . This reformulation has both conceptual and practical advantages; it expresses real- 
time response in terms of objects with intrinsic physical and geometric meaning. This is, for 
example, essential for the proof in [ji]] that various boundary transport coefficients can be 
solely expressed in terms of quantities evaluated at the horizon. 



In this note we further shed light on the formulation of [lOj] by showing that it follows 
simply from an analytic continuation of the original Euclidean formulation of AdS/CFT. At 
the end of our investigation, we came to realize that a similar procedure has in fact been 
used earlier [3] to justify the prescription of Q]. Nevertheless, we feel that our presentation 
may still be of some value, as it highlights the importance of the canonical momenta as the 
response of the dual operator even in real-time. 

We then show that our formalism may be extended to spinors, which involve interesting 
subtleties as we are now dealing with a first order system. We work out an explicit prescription 
for calculating real-time retarded correlators of fermionic operators, which is important, e.g., 
for probing quasi-particle structure associated with Fermi surfaces from gravity 111. 1 121]. 
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The plan of the paper is as follows. In section [TT] we state and derive our reformulation 
of the real-time prescription for retarded Green's functions in AdS/CFT. In section nil we 
work out in detail the application of this prescription to spinor fields in an asymptotic AdS 
spacetime. In section [IV] we further illustrate the prescription using two exactly solvable 
examples: pure AdS and BTZ black hole. We conclude in section |V] with a brief summary. In 
Appendix [Al we state our conventions for various Euclidean and retarded Green's functions, 
while Appendices |B] and O contain some technical details related to the main text. 



II. DERIVATION OF A REAL-TIME PRESCRIPTION FOR RETARDED COR- 
RELATORS IN ADS/CFT 

In this section we first state and then derive our reformulation of the real-time prescription 
for retarded Green's functions in AdS/CFT. For illustration, we consider a scalar operator 

which is dual to a massless bulk scalar field <fi. Massive modes contain extra divergences; 
while harmless, these complicate the discussion and are discussed in Appendix |B] The 
generalization to tensors is self-evident, and the generalization to spinors is discussed in 
detail in the next section. 

To be specific, we will consider a scalar action 

S =~\j d d+l x^-g{d4>? + --- (1) 
on a background spacetime metric of the form 

ds 2 = —g u dt 2 + g rT dr 2 + gudx 2 = g rr dr 2 + g^dx^dx" . (2) 

The boundary is taken at r = oo, where various components of the metric have the asymptotic 
behavior of AdS with unit radius, 

9tt,g rr ,gu ~ r 2 , r -> oo . (3) 

We will assume that the above metric has a horizon of some sort (to be more explicit below) 
in the interior. 1 

1 In a spacetime which does not have a horizon (e.g. global AdS), i.e. completely regular in the interior, (f> 
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We also assume that the theory is translationally invariant in x M directions (i.e, all metric 
components depend on r only), and work in momentum space along these directions, e.g. 



<f>(r, a?) = 0(r, k,) e --*+**, = (-co, k) . (4) 

We will frequently need to analytically continue to Euclidean signature via 

t — > —ir, uj — * iu)E, iS — > —Se (5) 
where the Euclidean action Se is given by 

Se = U dd+lx V9(d<P) 2 + --- ■ (6) 

A. Prescription 

We begin by simply stating the prescription of for computing the retarded correlator 
G^kfj) for O. In momentum space this reads: 2 

U(r:k„ 



G R {K) = lim 



j '"id 



(7) 

00=0 



<f>n(r, hp) 

where II is the canonical momentum conjugate to <fi with respect to a foliation in the re- 
direction. 4>R(r, k/j) is the solution to the equations of motion for </> which is in-falling at 
the horizon and satisfies the boundary condition lim^oo (j)R(r, kfj) — > 0o(^t)- The notation 
II(r; k^)\^ R in equation (j7|) indicates that II should be evaluated on the solution cj)R. Finally, 
the subscript 0o = means that in evaluating this ratio one should only take the part that 
is independent of <po', in an interacting bulk theory II and <pR will typically contain higher 
powers of <po which are not relevant for two-point functions. 



has a discrete spectrum. The corresponding boundary theory retarded function for O can then be written 
as a discrete sum of delta functions and thus does not have an interesting analytic structure. 
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The reason we have (J7]) instead of G/j(fc M ) = (lim 



is due to subtleties in taking the limit 

(f>o=0 r— 1 

r — > oo. An example is the massive scalar case discussed in Appendix [Bj Also, we found in [l0( that the 
fact that all quantities in ([7]) remain meaningful in the bulk at arbitrary radius is helpful for a physical 
interpretation of this formula. 
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Note that (j7|) can also be equivalently (in fact more generally) written as 

(0(M>0o= limn(r;^)U fl (8) 

r— >oo 

where ^{ky))^ denotes the response of the system to external perturbations generated by 
adding J d d x <p (x)O (x) to the boundary theory action. Equation (JTj) is the linearized limit 
of flE}. 3 

An alert reader will immediately recognize that is the exact Lorentzian analog of the 
standard prescription for boundary correlation functions in Euclidean signature as formulated 
m . Indeed we will now show that ([7]) and ([8]) (at the linearized level) can be derived 
from the original prescription of using a simple analytic continuation, even though they 
do not appear to follow directly from an action principle themselves. 



B. Analytic continuation 

First recall that a retarded correlator is analytic in the upper half complex-a; plane, and 
that the value of Gr(u, k) along the upper imaginary w-axis gives us the Euclidean correlator 
G E , i.e. 4 , 

G E (u E , k) = G R (iu E , k) u E > . (9) 

This expression applies both for finite and zero temperature, and both for bosonic and 
fermionic operators, with uj e only taking discrete values at finite temperature. Equation (Q 
can now be inverted to obtain Gr(uj, k) as 



Gr(uj, k) = G E (u E ,k) 



(10) 

u>E =— i(u}+ie) 



3 To sec this more explicitly, it is helpful to recall the standard result from linear response theory that if one 
considers a system in equilibrium at t — > — oo and then perturbs its action with the term J d d x 4>o(x)0(x), 
the one-point function of O in the presence of the source is given (to first order in 4>o) by 

(0(M}^o = GRikjMh)* 

where Gr is the retarded correlator of O. 

4 Ge{ue) for uje < is obtained from the advanced function. See Appendix \X\ for our conventions on the 
definitions of Euclidean and retarded functions. 
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Note that at finite temperature loe only takes discrete values, and so the analytic continuation 
can be tricky. We will simply ignore this issue below, assuming that Gr is sufficiently well- 
behaved that such an analytic continuation is possible. 

We now look at the central prescription of AdS/CFT in Euclidean signature (for 
both zero and finite temperature) 



exp 



d d x (f> (x)O(x) 



— g-<Sgravl<?><)J ^-Ql 

QFT 



where >Sgrav[0o] is the bulk action for <f> (which may include necessary boundary terms) eval- 
uated at the classical Euclidean solution <pE which is regular in the interior and satisfies 
boundary condition lim^oo 0£?(r; x) — > (f>o(x). From (TTTT) . one finds that one point function 
of O in the presence of the source <p can be written as 

<°to>*> = ~^rs = - lim n *(^)U (12) 

and its Fourier transform 

{O{uj E ,k))^ = - lim UEi^UEyk)^ . (13) 



In ( 1121) . is the canonical momentum for <fi in Euclidean signature and should be evaluated 
at the classical solution <pE- The second equality of (fl2"|) follows from the well known fact 
in classical mechanics that the derivative of an on-shell action with respect to the boundary 
value of a field is simply equal to the canonical momentum conjugate to the field, evaluated 
at the boundary. Note that in general the limit r — > oo in ( [TBI) is non-trivial and requires 
careful renormalization, which was developed systematically in {13I [jjj] . Similar cautionary 
remarks apply to other formulas below (in both Lorentzian and Euclidean signature). 

At linear level in <p , one then finds that 

U E (r, u E ,k)\^, 



G e (uJe, k) = - lim 



r— *oo 



(14) 

<po=0 



<j> E (r, uJ E ,k) 

In Euclidean signature, the above procedure is unambiguous as 0^ is unique (at least for suf- 
ficiently small O ); there will be only one solution that both satisfies the boundary conditions 
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at infinity and is regular in the interior. This is not the case in Lorentzian signature, where 
in the interior one typically finds two oscillatory solutions, both of which are regular. 

Once GE(^E,k) is obtained, one can of course obtain Gr using (flOl) . Here, we would 
like to argue that one can in fact obtain an intrinsic prescription for retarded functions 
by analytically continuing (f>E{r,WE,k) directly in the bulk to Lorentzian signature 5 . More 
explicitly, introducing 

4>r(uj, k) = 4> e (uje, k) 



(15) 



and then from (TTOT) and (fMl) . we obtain simply' 

n i 1 \ t ^-E{r,uJE,k)\^ 
G r {uj, kfi) = - hm 



r— >oo 



= hm n(r; "' fc) K (16) 
u E =-i{u+ie) r ^°° (p R {r,uj,k) 



<t> E {r, u E , k) 

which gives (J7J). Similarly, ( TT3T) becomes (jHJ) under the same analytic continuation. To 
complete the derivation we still need to show that 4>r obtained in (ITBT) satisfies the boundary 
conditions at the horizon and infinity as stated below (J7j). Since analytic continuation does 
not change the boundary conditions at infinity, we need only show that the function on the 
right-hand side of (ITS"]) satisfies the in-falling condition at the horizon. We do this explicitly 
by examining the behavior of <pE near different types of horizons 7 : 

• N on- degenerate horizon: consider a Euclidean black hole background of the form 

ds 2 = f(r)dr 2 + -r^dr 2 + a{r) 2 dx 2 (17) 
f{r) 

where f(r) has a simple zero at r = r , i.e. f(r) ~ -w(t — ro) and r is periodic in (3. 
From the wave equation of 0, one finds the following near-horizon solutions 

4>+(r) ~ (r — ro) 4?r or 4 > -{ r ) ~ ( r — r o) **~ ■ (18) 
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A similar analytic continuation has been used earlier in to justify the prescription of 

6 Note the relative minus sign between the first and second equality below is due to that in Euclidean 
signature Ii£ = y/gg rr d r cf> (see ([6])) and in Lorentzian signature II = —^—gg rr d r (j) (see ([1])) with an extra 
minus sign. 

7 Note that the following discussion applies to any mass, as the mass term is never important at the horizon 
with a nonvanishing oj. 
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When the real part of uje > (which is always the case when evaluating a retarded 
correlator by (TTOil ). ~ (f> + and thus we find that 



~ (r — r ) ^ (19) 

oj£ = — i(oj+ie) 



which is precisely the behavior of an in-falling wave, since from (J3J) it leads in coordinate 
space to a wave of the form e _lw *(r — tq)~^ . 

Degenerate horizon: For a degenerate horizon (e.g. of an extremal black hole), we 
have a metric of the form ( TlTl) except that / now has a double pole at r = ro, i.e. 
/ ~ c(r — r ) 2 , where c > 0. The discussion is similar to that above except that the 
near-horizon solution can now be shown to take the form 



4> e (uje, k) 

again corresponding to an in-falling wave 



~e c ('— -o), (20) 



Poincare horizon: finally let us consider a Poincare horizon of the form (near z — > oo) 

ds 2 = — \dr 2 + dx 2 + dz 2 } (21) 
z 2 

Then the near horizon solutions for <fi are given by 



<f>± ~ e ±kz , k = yj u% + fc 2 , z -> +oo (22) 

Note that if we take the branch of the square root that is positive, then the regular 
solution <j) E ~ 4>-. We now perform the analytic continuation u E = —i(uj + ie), under 
which the branch of the square root with positive real part becomes 



iy/ uj 2 — k 2 uj > \k\ 



(23) 



iy uj 2 — k 2 uj < — \k\ 



where we have assumed a timelike momentum \uj\ > \k\ (as otherwise the problem is 
the same as in the Euclidean case). When plugged into 0_ ~ e~ kz , Equation (|23|) again 
describes an in-falling wave for both signs of uj. 
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The key point here is that the +ie in the analytic continuation (1101) is important; it guarantees 
that the solution that was regular at the horizon in Euclidean signature becomes the solution 
that is in-falling at the horizon in Lorentzian signature. 

Note that while the Euclidean prescription (|TT]) - (|T4|) for computing correlation functions 
follow from an action principle, this is not the case for the Lorentzian prescription ([7]) and ([8]), 
which cannot be obtained by taking functional derivatives of the on-shell Lorentzian action. 
While analytic continuation of a Euclidean solution yields a Lorentzian solution as in ( TT5T) . 
analytic continuation of a Euclidean on-shell action-which necessarily involves an integral 
over the full spacetime manifold-does not necessarily yield the correct Lorentzian action. The 
Lorentzian manifold generically differs from the Euclidean one in many crucial aspects such as 
the topology, the number of boundaries, etc., making the implementation and interpretation 
of such a procedure much more subtle. On the other hand, the analytic continuation of the 
canonical momentum as in our prescription requires only a local analytic continuation on the 
boundary of the manifold and is thus much simpler. 

C. Vector and Tensor Operators 

Before concluding our discussion of bosonic operators, we point out that in this formalism 
the tensor structure of the real-time response follows naturally without needing to decompose 
the system mode by mode into a series of scalar wave equations. Let us imagine first a bulk 
vector field Am with Maxwell action and gauge coupling g e fr; this is dual to a conserved 
current J' 11 , and the expression ((HI) then becomes 



where the expression on the right-hand side is the momentum conjugate to A^. Similarly, if 
we consider the field theory stress tensor T Miy , we find that 



correctly captures the real-time response of the system to metric perturbations. Here 7^ 
and Kn V are the induced metric and extrinsic curvature on each constant-r slice, and the 




(24) 




N 



(25) 



9 



161 ]. i.e. precisely the 



expression on the right-hand side is the Brown- York stress tensor [15 
momentum conjugate to a gravitational perturbation h^. The validity of these formulas even 
in real-time was a key element in the results relating AdS / CFT to the membrane paradigm 
in M. 



III. RETARDED CORRELATORS FOR FERMIONIC OPERATORS 

In this section we generalize the prescription (]7j) and (jHJ) to fermionic operators. We first 
consider the calculation of Euclidean functions in an asymptotic AdS geometry, reviewing 
(and slightly generalizing) the earlier results of {17I . 1^ . 19]. We then discuss the prescription 



for real-time retarded correlators. In the next section we use the prescription to calculate 
retarded correlation functions for fermonic operators in field theories dual to the pure AdS 
and BTZ geometries, where closed-form expressions can be obtained. 

We consider a boundary theory fermionic operator O which is dual to a spinor field tp in 
the bulk 8 . We will suppress all spinor indices. Since we are interested in two point functions 
of O, it is enough to consider the quadratic part of the action for ip, which can be written as 

S =Af J d d+1 x^gi{i[)T M D M iP - rmfjip) + S hd (26) 

where ip = ip^T* and 

D M = d M + l -u abM T ab . (27) 

denotes the boundary terms required to ensure that the total action has a well defined 
variational principle [3] and is briefly discussed in Appendix O In (I26j) - (127j) M, N ■ ■ ■ 
denote abstract spacetime indices and a,b, - ■ ■ denote abstract tangent space indices. Gamma 
matrices with a specific index (like r*, T r , P as opposed to those with an abstract index T M ) 
always correspond to those in the tangent frame. The boundary theory gamma matrices will 
be denoted by 7 M . The background metric is taken to be (j2J) with the associated asymptotics. 



Note that in this section O is fermionic whereas in Section [IT] O was bosonic; we apologize for the abuse of 
notation. 
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We will be working in momentum space with the Fourier transform of tp denoted by if)(r, k^). 
Finally, the action is normalized by a factor A/"; we will ignore this factor in the following, 
as it simply contributes an overall factor in front of all boundary theory correlators, but one 
can show that its sign is fixed by bulk unitarity. 9 

The analytic continuation to Euclidean signature is as in (J3J) with also 

7* -> -iY, T* -> -iT T , i) -> (28) 
and thus the Euclidean action following from (|26p 

S E = ~ [ d d+1 Xy/g{i>Y M V M ^ - m^) + S hd . (29) 



A. Euclidean correlators 



We first review how the Euclidean prescription flli p — (j!4p works for a spinor operator, 
where we now have 



cxp 



d d x {xoO + Oxo) 



— g-' S 'grav[X0 1 X0] (30) 

QFT 



To find the right hand side of (130|) . we need to construct a ipE which is regular in the 
interior and satisfies the boundary condition lim^oo ipE = Xo- Attempting to interpret this 
equality, we realize that the situation here is a bit more subtle than for a scalar field: ip 
and Xo are spinors of different spacetime dimensions, and thus may have a different number 
of components. Also, the action for ip contains only one derivative and imposing Dirichlet 
boundary conditions for such a first order system requires more care. It turns out that these 
two issues are intimately related. 



9 To see this, consider canonically quantizing the action (|26[) by imposing the equal-time anticommutation 
relation {ip(x), n(x')} = i5^ d \x, x'), where here n is the momentum conjugate to ip with respect to time. We 
find that 7r = —iNy/—gip^. If we now require that the anticommutator {x/j, tp'} be positive~& requirement 
in a Hilbert space with states of positive norm-then we find that this fixes the sign of Af to be negative. 
Reinstating this factor of Af in front of retarded correlators found in this paper results in positive spectral 
densities in the boundary theory, as required by boundary unitarity. 
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To see this, it is convenient to decompose ip in terms of eigenvalues of T r , with 

^ = ^ + + V-, ^± = r ± ^, r± = ^(i±r) (31) 

where from (1291) we find the corresponding Euclidean canonical momentum 10 IT-|- (in r-slicing) 
conjugate to ip± to be 

_i _ _ _i _ 

n+ = -y/ggrr 2 ^-, n_ = ^gg rr 2 if) + . (32) 

Thus we see that ip± are conjugate to each other. If we supplied Dirichlet boundary conditions 
for both of them at infinity, then it would completely fix the solution everywhere; however 
this is obviously incorrect, as (generically) this solution would not be regular in the interior. 
Instead we should begin by demanding regularity in the interior; this then leaves us with 
the freedom to impose boundary conditions Xo f° r either or ip_, i.e. only for half of the 
components of ip. 

When the boundary theory dimension d is even, we can choose the bulk Gamma matrices 

= y, r r = Y +1 (33) 

where >y d+1 is the analogue of 7 5 for d = 4. Thus from the (/-dimensional point of view, the 
two components ip± transform like d-dimensional Weyl spinors of opposite chirality. As the 
boundary value of one of ?p±, Xo (and so also O) is then a (/-dimensional boundary spinor of 
definite chirality. Thus for d even a Dirac spinor ijj in the bulk is mapped to a chiral spinor 
O on the boundary. 

When d is odd, it is convenient to work with the representation 



r = 



1 > 




'07M 




1 r H 


\o-i) 




1 r J 



(34) 



which can be easily seen to satisfy the d+1 dimensional Clifford algebra. In this basis we 
can see that the two components tjj± each transform as a d-dimensional Dirac spinor; thus 

10 For fermions we use only the Euclidean canonical momenta, and have thus omitted the subscript E 
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for d odd Xo an d O are both Dirac spinors. In all dimensions the number of components of 
O is always half of that for ip. 

To decide on which of i/j± to impose Dirichlet boundary conditions, we must examine their 
asymptotic behavior at large r, which can be worked out by solving the Dirac equation for 
ip in the region r — > oo (with metric given by (j3J)). This computation is performed in more 
detail in Section TlV Al and here we discuss only the asymptotic behavior of i/j, which is given 
by 

ip+(r,k) = A(k)r-% +rn + B(k)r-%- m -\ </>_ = C^r"^™" 1 + D(k)r~^- m , r -> oo 

(35) 

Plugging this expansion back into the Dirac equation we find the following relations between 
the expansion coefficients 

D = -^(2m+l)B, C=^J^A, 7 • A; = 7%, k 2 = k^ . (36) 

Using ( l32i) . we see that the corresponding canonical momenta behave as 

11+ = -Cr^™- 1 - Dr^- m , n_ = Ari +m + Br^ m ~\ r -> 00 . (37) 

Note that in deriving (|36|) we have used (1331 or (|34|) for d even or odd. 

Notice that as we take m — > —m, we simply exchange the role of ip±, with A ^ D and 
B ^ C. We can thus restrict our attention to m > 0. Among all terms in (|35l) . the term 
with coefficient A is dominant. We thus should the impose boundary conditions 

A = xo, i-e. lim A~ m ^ + = X o ■ (38) 

r— >oo 

Now taking a derivative with respect to xo on both sides of (1301) . one finds that as in ([12]) . the 
response of O is formally given by the conjugate momentum n + . n However, for general m, 
as in the case of a massive scalar discussed in Appendix (Bj one should define the boundary 
limit carefully. In view of the second equation in ( l38l) . we thus find 

(O) xo = - lim r m -fn + , i.e. (O) X0 = D. (39) 

r— »oo 

11 See Appendix [Cl for an explicit discussion of the variation of iSgrav 
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Note that in obtaining the second equation above we have extracted only the finite terms in 
the right hand side of the first equation 12 . This appears reasonable since the other term in 
which is proportional to C, is locally related to the source A. 

With this boundary condition on ip(r — > oo) specified, the solution iJje to the equations of 
motion which is regular in the interior is then uniquely determined. From this we can extract 
the corresponding D; we will find that D and A are related by a matrix S, 

D(k) = S{k)A(k) (40) 

then the boundary Euclidean two-point correlator is given by 13 

G E (k,) = S(k,)Y . (43) 

Sometimes it may be easier to solve for B in fl35l) . in which case introducing a matrix T 

B = TA (44) 

we then have from (136j) 

G E (k ll ) = -^(2m+l)( 1 .k)TY ■ (45) 
To conclude this subsection, we make some further remarks: 



1. Using the standard scaling argument, the identification of the source and response 
in (I38p and (|39|) implies that the scaling dimension A of O is related to m by 

A = ^ + m (46) 



which is consistent with results obtained in 



171 ] for pure AdS. 



12 Here we are assuming that divergent terms should be removed by holographic renormalization. It might 
be worth checking this more explicitly. 

13 Note that from 

(0(x)) = J d d yG E (x~y)^ Xo (y) (41) 
which in momentum space becomes 

(0(k)) =G B (fc)7 T Xo(fc) (42) 
The 7 r factor is due to that G E ~ (OO t ) rather than (OO). 
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2. In the case of a massive scalar (as discussed in Appendix [B]), the r — > oo limit in 
the ratio -P- is somewhat subtle since subdominant terms in II and &e also contribute 

<t>E Tr> 

which changes the overall constant. This does not appear to happen to fermions. 

3. When < m < |, all terms in ip± are normalizable. We thus can choose either A or D 
as the source and treat the other as the corresponding response. Note that if we choose 
D as the source term then 

d d — l.d , lf , s 

A = m, ->■ < A < - (47) 

2 ' 2 2 v ; 

In this range the double trace operator of O is a relevant perturbation. The partition 
functions for the two alternative ways to quantize the theory should be related by a 
Legendre transform as one is the conjugate momentum of the other j^ . 

4. For m = |, the two terms in are degenerate. Instead one has 

ifj_ = r ~2-\(C log r + D) (48) 

The relations between A, C in (|36|) are now replaced by 

C=(iy k)A . (49) 

In this case the term proportional to A is not normalizable and should be treated as 
the source term. 

5. For even d, if Xo has negative chirality, then we should be imposing boundary conditions 
on if)_, which in turn implies that the bulk mass term must be negative. An example 
illustrating this is shown in Section TlV Bl with the BTZ black hole. 



B. Prescription for retarded correlators 

Now we simply analytically continue ipE and (|38|) - (|45|) to Lorentzian signature, precisely 
as in the bosonic case with ffl5l) and (flOl) . We obtain the following prescription for computing 
retarded correlators for spinor operators: 
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1. Find a solution ^(r, A;^) to the Lorentzian equations of motion which is in-falling at 
the horizon. 

2. Expand ipFi(r, k^) near r — > oo as 

ifj R+ = Ar~i +m + Br"-*-™- 1 , rjj R _ = Cr"^" 1 " 1 + Dr"^ m , r -> oo (50) 
where ^± are spinors of definite r r eigenvalue as defined in (j3TT) . 

3. Then Gr^) is obtained from the analytic continuation of (jl3"]l by 

9rr) -I- 1 

= <5(fc M ) 7 « = • *) (51) 

where 7 • k — 7^^, A; 2 = fc^A;^ and 5 and T are defined by the relations 

D = SA, B = TA . (52) 

We end this section by noting that causality and unitarity of the boundary CFT impose 
two important properties on its real-time correlators. Causality implies that the retarded 
correlator must be analytic in the upper half of the complex uj plane; this is easily seen 
to be true for correlators calculated using the prescription above. Unitarity implies that 
the imaginary part of the diagonal elements of the matrix Gr is proportional to a spectral 
density and must be positive for all u; as mentioned above, the overall sign of the correlator is 
proportional to the normalization of the Dirac action, and one can show that a sign consistent 
with bulk unitarity also results in positive spectral densities in the boundary theory. 



IV. EXAMPLES 



In this section we apply the prescription developed in last section to two exactly solv- 
able examples, fermionic retarded correlators from pure AdS and a BTZ black hole back- 
ground [2 11 ]. 



16 



A. Pure AdS 



The Euclidean vacuum two-point function for a spinor operator O in a CFT, which is dual 
a fermonic field ip in a pure AdS, was obtained in closed form before in 17J] (see also 18l. Il9|). 
The corresponding retarded function can then obtained from it by analytic continuation 
using ({TO]) . Here we calculate the retarded function directly using the prescription developed 
earlier check of our formalism. 



For pure AdS, the metric is given by 



ds 2 



-dt 2 + dx 2 ) + 



dr 2 



+00 . 



with associated nonzero spin connection components given by 



Utr 



—rdt, Mi, 



rdx 1 



The Dirac equation in momentum space is then given by 

rT r d r ip + -T ■ kip + -T r ip - rmp = 
r 2 

Using ( 13TT) . equation (1551) now becomes coupled equations for i/)± 



0H 



i-y ■ k 
k 2 



A(—m)ip_ 



0- 



—g r A(m)ip+ 



where we have used ( 133]) or ([Mil for d even or odd and 



A(m) = r I r<9 r + 



d 



m 



from which we obtain 



k 2 ip + = A(—m)A(m)ip + . 



(53) 



(54) 



(55) 



(56) 



(57) 



(5* 



Note that ( |58l) is now a group of decoupled scalar equations, which implies that the T matrix 
defined in fl52|) will be proportional to the identity matrix. Equation fl58|) can be solved exactly 
using Bessel functions, and the solution satisfying the in-falling boundary condition at the 
horizon is 

( / \ 

k 2 >0 



r 2 ^m+i 



r 2 



™ 1 1 I r- 



_±ti (2) 

m +2 



^UJ 2 -\k\ 2 



a + a> > |fc| 



a + w < — Ifcl 



(59) 
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where a + is an arbitrary constant spinor. The story for the spacelike case k 2 > is exactly 
the same as that of Euclidean correlator and the retarded correlator is real. For the timelike 
case u > \k\, we find the corresponding A, B coefficients as defined in (150}) as 

giving 

V = vfef(^) 2 " ,+1 e- < ™ + ^ <«) 



Using ( 1511) . we then find that 



g fl (ib)= ' ' 1 ^ (7^)^ *>|£| (62) 



k 2 r(m + |) 

For the other possible timelike case ui < —\k\ we simply change the phase factor e~^ m+ ^ 7 

to e ( m +!K 



B. BTZ Black Hole 



We now consider fermionic correlators in a BTZ black hole background 21|. Previous 



work 



221 ] found that the quasinormal modes of a BTZ black hole precisely coincide with the 



theory 



22 



(whose form 



23] to solve the 



poles in the retarded propagator of the appropriate operator of the dual 
is essentially fixed by conformal invariance). Here we will closely follow 
wave equation but slightly extend these results by finding the full correlator from the gravity 
side. 



A BTZ black hole with a mass M and angular momentum J describes a boundary 2d 
CFT in a sector with 

L = l -{M + J), U = \{M-J). (63) 

The system has a finite entropy and non-vanishing left and right temperatures. Here we are 
interested in probing this sector using fermionic operators with spin s = \. Such an operator 
0± is characterized by conformal weights (hi, Hr) with 

1 



h T , + h 



R 



X 



hi, - hi 



±- 



(64) 
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where the ± sign denotes its chirality. As described in the previous section each + (or 
O ) is described by a Dirac spinor if) in the bulk, with different chiralities corresponding 
to different boundary conditions and opposite bulk mass. We will now find the retarded 
correlators of 0± by solving the Dirac equation for if) in the BTZ geometry. 

The metric of a BTZ black hole can be written as 

(r 2 - r 2 )(r 2 - r 2 ) 2 r 2 dr 2 2 / r + r_ \2 

ds = ~2 dt + ( r 2_ r 2 )(r2 _ r 2 } + r {^-—M) ( 65 ) 

where <fi is an angular coordinate of period 2n. The mass, angular momentum, and left and 
right moving temperature of the system are given by 

M = ^G^' J= ^' Tl = —^T Tr = —^T (66) 

where G is the 3d Newton constant. To solve the Dirac equation in (168|) . it is convenient to 
switch to a new coordinate system (p, T, X) 

r 2 = r\ cosh 2 p-rl sinh 2 p, T + X = (r + + r_)(t + 0), T-X = (r + -r_)(t- 0) (67) 

in which the metric is 

ds 2 = - sinh 2 pdT 2 + cosh 2 pdX 2 + dp 2 (68) 
and the spin connections are given by 

u)Tp = — cosh pdT uxp = sinh pdX . (69) 

We will work in Fourier space on each constant-p slice; a plane wave can be decomposed 
in either the (X, T) or (0, t) coordinate system: 

tp = e- lkTT+lkxX ifj( Pl A; M ) = e- Mk<t >i)(p, fc„) (70) 

where using (1671) and (|66|) we can see that the momenta in the original coordinate system 
(u, k) are related to (fc^, ^x) by 

j j u + k uj — k 

k T + k x = ——, k T -k x = ——. (71) 
2irl R Itx1 l 
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The Dirac equation can then be written as 



f 1 /cosh p t sinh p 



k x T x 



k T T T 



— m 



2 ^sinhp coshpJJ \yCoshp sinh p 
We choose a gamma matrix representation where T p = a 3 ,T T = ia 2 ,T x 



if) T = (ip + ,ip-). Now following 23[ and letting 



cosh p ± sinh p 



cosh p sinh p 

then in terms of Xi,2 and 2, the Dirac equation becomes 



z = tanh 2 p , 



V? = (72) 
a 1 and write 

(73) 



2(1 - z)yGd z xi - i \Jj= + kx^/z^j xi= (rn - ^ + i(fcr + k x )*j 



^ + i{k T + k x ) ) X2 

2(1 - z)Vzd z X2 + i[-$= + kxVz] X2 = ( m - \ - i(k T + kx) ) Xi 



(74) 



Note that the horizon is at z = and the boundary at z = 1. It is now possible to eliminate 
one of the fields Xi,2 to obtain a second-order equation in the other field. The solutions to 
that equation are given in terms of hypergeometric functions. We are interested in evaluating 
the retarded correlator, and so we pick the solutions that are infalling at the horizon. These 
take the form 14 



X2(z) = z a (l-zfF(a,b; c; z) 
,' a — c\ 1 

Xi{z) 



zi +a {l- z) F{a,b+l;c+l;z) 



where the parameters are 



ikx n 1 m 

a = 3= 1 

2 H 4 2 



and 



a = 2 ( m+ \) ~ 2~^ T ~ kx ^ b= \\ m ~\)~^ kT + kx ^ c= \~ ikT - 



(75) 



(76) 



(77) 



14 In order to display consistency of these solutions with the equations of motion ([74| . it can be helpful to 
use the hypergeometric identity —aF(a + 1, b + 1, c + 1, z) + j-~^F(a, b, c, z) + j^-F(a, b + 1, c+ 1, z) = 



24] 
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From (173 p . we know that at the AdS boundary tp has the following asymptotic behavior 

A(l - z)~- + B{\ - z) 1+ ~ ~ C(l - z) l ~- + £>(1 - , (78) 

Our earlier analysis tells us that if m > we can identify A = xo as the source and .D = 
as the response, and so the retarded correlator in this frame is given by 

Gr = (79) 

We now explicitly expand the solutions ( 1751) near the boundary to extract the coefficients D, 
A. The relevant ratio works out to be 

r (k k ) - F (j ~ ™) r (K 1 ~ 2 ^ ~ ^) + 2m )) r (|(3 - 2z(fc T + fo) + 2m)) 
R[ T} X) ~ V (| + m) r (1(1 - 2t(fcr + A*) - 2m)) T (1(3 - 2t(k T - k x ) - 2m)) 

(80) 

Using the relations flTTj) this can be rewritten in terms of momenta in the (t, <fi) coordinate 
system as 

g a = -i r ); -) £ ( ) m { (si) 

r( 3 + »n)r(^-i^)r(^-i^) 

where we have introduced 

, m 1 , m 3 . n . 

h L = h -, /ir = h - (82) 

2 4' 24 v 7 

and 

~ m 3 ~ m 1 . 

^ = -2+4. ^ = "Y + 4- (83) 

Note that ( 1HT1) has a nice factorized form for left and right sectors. The correlator has a pole 

whenever either of the two gamma functions in the numerator has an argument that is a 

negative integer; thus we find the following two sequences of poles 

u — -k- A-KiT R (n + h R ) u = k - 4iriT L (n + h L ) n e Z + (84) 

As demonstrated in 



22J, these two sequences of poles are precisely those appearing in the 
finite-temperature retarded correlator of an operator in a 2D CFT with conformal weights 
(hL, hii). Since — Hr = — |, this is consistent with our expectation that tp corresponds to 
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If m < 0, we then find that tft_ is the source and ip + is the response. Thus we are looking 
the correlation function of 0+. One can immediately find the correlator from the results 
above, as the roles of D and A are reversed; we find that this results in a different pole 
structure consistent with the correlator of an operator of conformal dimensions (h^, /i#), 
indeed corresponding to a boundary spinor of positive chirality. 15 

Finally we note one last point which is important if one wishes to determine the overall 
normalization of the correlator. Note that the above expression was computed in the (X, T) 
coordinate system; however we see from (IBTj) that the relation between this and the (t, (ft) 
coordinate system involves scaling the left-moving coordinate by (r+ — r_) and the right- 
moving coordinate by (r+ + r_), i.e. by the left and right temperatures respectively. Thus 
the correlator in the (t, (ft) frame is 16 

G R (u, k) = (27rT L ) 2h ^ 1 (27TT R ) 2h ^- 1 G R (k T , k x ) . (85) 

This rescaling is important for the extremal limit — > 0. In this limit, using the Stirling 
formula in ([8*T1) . we find that G R blows up 

r (I + m ) r (h R - ig±|) V l J 

but G/j does have a finite limit 



r (I + m) r ^ K _ 



— ^ 



'4ttTr y 

This expression can be verified by direct calculation using the bulk spacetime for the extremal 
BTZ black hole with = 0. Note that (!H7|) again has a nice factorized form with the left- 
moving sector given by the vacuum expression. 



15 



Note that in attempting to compare our assignment of conformal dimensions directly with [221 ] , one should 



keep in mind two issues: m in this paper is —m in [221 ] . and our choice of gamma matrices means that 



a 3d spinor ip + with positive eigenvalue of T p has opposite 2d hclicity here than it does in 22| ; thus our 
assignment of Hl and is switched relative to them. 
16 — l's in the exponents of the expression below are due to the fact that this is an expression in momentum 
space. 
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V. CONCLUSION 



In this paper we showed that an intrinsic Lorentzian prescription for computing retarded 
Green functions from gravity can be obtained by analytic continuation from the corresponding 
problem in Euclidean signature. An important message here and in our earlier work 10] 
is that, even in Lorentzian signature, for any bulk field- whether bosonic or fermionic-the 
corresponding conjugate momentum contains the response of the dual operator. The field 
theory response can thus be expressed in terms of quantities with clear physical meaning 
in the bulk, making more transparent the origin of phenomena such as the universality of 
transport coefficients [lo| . 

We also explained in detail the issues that arise when attempting to find Lorentzian cor- 
relators of fermionic operators, and worked out simple examples to demonstrate the method. 
It is easy to think of future applications of these methods. In particular, there has been 
relatively little investigation of fermionic correlators at finite temperature or density using 
gauge-gravity duality (for recent work see 25] ), and it would be interesting to see whether 
effects such as the universality of transport coefficients have fermionic analogues. Recent 
work with the real-time response of spinor operators [3] has also found new phenomena 
such as the existence of a Fermi surface at finite chemical potential, and it remains to be 
seen whether fermionic probes can help reveal other unexplored structure in gauge-gravity 
duality. 
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APPENDIX A: DEFINITIONS OF CORRELATORS 



We start by assuming that O is a hermitian bosonic operator; in that case our convention 
for the Euclidean correlator is 

G E (t,x) = {T E O(x)O(0)), (Al) 

where T E denotes Euclidean time ordering and at finite temperature the Euclidean time 
direction is taken to have period (3. 

Our conventions for the various realtime correlators are as follows, where p denotes the 
thermal density matrix: 

G R (t,x) = ie(t)tr(p[O(t,x),O(0)}) (A2) 

G A (t,x) = -ie(-t)tr(p[O(t,x),O(0)}) (A3) 

Note that with this sign convention the imaginary part of Gr is equal to n multiplied by 
the spectral density, and thus is positive definite. Another common convention is to use the 
opposite sign for both Gr and Ga\ in particular, this was used in our earlier work 



For O fermionic and complex we use the following conventions 

G E (t,x) = (T E O(x)O\0)), (A4) 

G R (t, x) = ie(t)tr(p{0(t, x), 0t(O)}) (A5) 

G A (t, x) = -i0(-t)tr(p{O(t, x), O\0)}) (A6) 



APPENDIX B: MASSIVE FIELDS 

For massive scalar fields, the boundary limit is more subtle due to various divergences. 
To be specific we consider a scalar action of the form 

S = ~\J d d+1 x^ {{3(f)) 2 + m 2 2 ) (Bl) 

with the background metric satisfying the standard asymptotic AdS behavior ([3]) near the 
boundary. 
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Then 0^ has the asymptotic behavior 

Mr, K) « A iK) rA ~ d + B ( k ») r ~ A i r ~> 00 (B2) 

where 

A = ^ + !/, z/ = ^™ 2 + j (B3) 

Note now the field theory source should now be taken to be A, which differs by a power of r 
from the boundary value of <p. This implies that 



{0(k,)) A = \]mr^ d U(r;k,)\^ (B4) 

r— >oo 



where as before II = —y/—gg rr d r (p is the canonical momentum to <ft. Plugging in the near- 
boundary expansion (1B2I) . we see that II has the asymptotic behavior 

n M/0L « —(A — d)A(k l j,)r A + AB(k fl )r d ^ A J (B5) 

and thus the extra power of r in (1B4I) exactly cancels the power of r in front of the sub- 
leading solution B, so the final answer contains a finite part that is independent of r. Note 



261 ] . and the retarded Green function 



nevertheless that some care must be taken in the limit 
Gr should now be written as 

G R {k,) = lim r 2 ^> n(r ;^« = (2A - (B6) 
r ~*°° (pR{r, ky) A^ky) 

where one is instructed to extract only the finite piece on the right-hand side. 



APPENDIX C: BOUNDARY TERMS FOR SPINORS 

Here we consider various technical issues related to the variation of the Dirac action. 
In particular, we explain why the momentum conjugate to can truly be considered the 
response of the operator (O) dual to ip; essentially this argument is due to [3] and we review 
it here for clarity. Throughout this section we assume m > 0; a similar argument holds for 
m < with an interchange of i/j±. 
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Consider the bulk Euclidean Dirac action 

Sbuik = - J d d+1 x^^T M V M ifj - mW) (CI) 
From this action the momentum Il + conjugate to ip + is 

n + = -Vgg^t- (C2) 

It appears that the momentum conjugate to ip is identically zero; however we expect this 
momentum to contain information regarding the field theory operator O, which certainly 
does not vanish in general. Upon an integration by parts, however, we can transfer the radial 
derivative from ifj to ■0; in that case we find that the momentum conjugate to ^+ is 0, but 
the momentum IT + conjugate to is now 

n + = -y=^V- (C3) 

Confusingly, however, it appears that the two expressions above cannot apply simultaneously, 
although our expectation from the field theory side is that neither of them should vanish. 
Clearly boundary terms in the action are playing an important role; we now explain precisely 
how to fix these boundary terms and what their effects are. 

The portion of the action (ICip containing radial derivatives can be written explicitly as 

Sbuik D-J d d+1 x^^(i;_d r ij + - iJ + d r ij.) (C4) 

From this we can see that if we vary this action around a solution to the equations of motion 
ip — > ip + dip; the variation works out to be 

SS bu0i = bulk term - / d d Xy/g^ FF {ip^5i) + - ^+^_), (C5) 

JdM 

where the bulk term is proportional to the equations of motion and the boundary term follows 
from flC4j) . Note that the on-shell action is thus a function of both ip + and ?/>_; however, this 
is not correct. As mentioned in the text, if we impose infalling boundary conditions we are 
no longer free to choose ?/>_, and this information must be implemented in our variational 
principle. This is done by adding to our action the following boundary term 



Sa = ~ / d d x^W^- (C6) 

'dM 

26 



We now see that the variation of the full action is 



5S 



total 



£(<Sbulk + Sq) 



(C7) 



dM 



Thus the on-shell action no longer depends on also, note that we have 17 

$<Stotal / — 7 T=r OiStotal 



n 



ih 



(C8) 



where we are now defining the conjugate momenta IT to be the on-shell variation of Stotai- 
Thus we see that the naive relations (IC2I) and (1031) do indeed hold simultaneously when we 
use the correct boundary action. It is easy to see that if we had started with a different 
action related to (IClj) by only boundary terms we would have obtained a different value of 
Sq but the final answer (IC8|) would have been the same. (It is shown in 19| that a symmetric 
splitting of the kinetic term in (|C1D results in the boundary term used in the original work 
with spinors in AdS/CFT [17j,[l8|]). 
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